Dynamics of Coupled Qubits Interacting with an Off-Resonant Cavity 
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We study a model for a pair of qubits which interact with a single off-resonant cavity mode and, 
in addition, exhibit a direct inter-qubit coupling. Possible realizations for such a system include 
coupled superconducting qubits in a line resonator as well as exciton states or electron spin states 
of quantum dots in a cavity. The emergent dynamical phenomena are strongly dependent on the 
relative energy scales of the inter-qubit coupling strength, the coupling strength between qubits and 
cavity mode, and the cavity mode detuning. We show that the cavity mode dispersion enables a 
measurement of the state of the coupled-qubit system in the perturbative regime. We discuss the 
effect of the direct inter-qubit interaction on a cavity-mediated two-qubit gate. Further, we show 
that for asymmetric coupling of the two qubits to the cavity, the direct inter-qubit coupling can be 
controlled optically via the ac Stark effect. 

PACS numbers: 78.67.Hc,85.65.+h,85.75.-d 



I. INTRODUCTION 

The quantum- mechanical coupling of individual atoms 
and photons has been demonstrated in a scries of cavity- 
QED experiments during the past decadei^In addition 
to the fundamental interest in states formed by super- 
positions of matter and photon degrees of freedom, the 
development of cavity-QED has been triggered by possi- 
ble applications in quantum information schemes. Strong 
coupling of a qubit and a cavity mode allows one to con- 
vert localized qubits into flying qubits suitable for trans- 
mission of quantum information^ Transmission through 
a detuned cavity mode can be utilized for qubit read- 
out For several qubits interacting with a single cavity 
mode, the cavity can act as bus for quantum information 
that couples qubits with large spatial separation. 5,7 

In this manuscript, we consider a pair of qubits both 
of which are coupled to one cavity mode and, in turn, 
both qubits are coupled directly to one another by an 
exchange-type interaction. The Hamiltonian of the sys- 
tem under study is given by 

H = [QiSi,z+9i(asf + aUr)] (1) 

T < - ■ +x 
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Here, Mli is the energy splitting of qubit i = 1,2, hgi 
is the coupling strength between qubit i and the cavity 
mode, Hu> is the energy of a cavity photon, sf are raising 
and lowering operators for qubit i, and Sj iZ = {\/2)a i z 
are pseudo-spin operators with the Pauli matrix a^ z . Ex- 
tending previous work 5 ^ on cavity-QED with solid-state 
systems, we take into account an additional direct qubit 
coupling with amplitudes HJj_, hJ z . The interplay of di- 
rect and cavity-mediated coupling leads to intricate phe- 
nomena. For J z = Jj_ =0, the state of both qubits can 
be read out simultaneously from the cavity dispersion if 
f2i 7^ f^2 and the cavity Q-factor is sufficiently largest 



We show that a similar readout scheme remains effective 
even for the case of coupled qubits. For identical qubits, 
fli = and g\ = g 2 , the measurement via the cavity 
mode projects the coupled state onto the singlet-triplet 
basis and thus would allow one to detect entangled two- 
qubit states. We show that the direct exchange coupling 
modifies the cavity-mediated interactions and decreases 
the fidelity of cavity-mediated two-qubit gates* 5 -' 8 On a 
more fundamental level, our model allows us to investi- 
gate the transition from cavity-mediated inter-qubit in- 
teractions at small J± to a system of two directly coupled 
qubits interacting with a cavity. 

An implementation of the Hamiltonian Eq. (TjQ) is fea- 
sible with various physical systems. For Jj_ = J z = 0, 
Eq. |TJQ) has been derived for superconducting qubits at 
the charge degeneracy point interacting with a microwave 
resonator £ In the present work, this scheme is extended 
to capacitively coupled superconducting qubits. Other 
implementations include exciton states in coupled quan- 
tum dots or coupled quantum shells in optical microcavi- 
tiesi&iS We show that this model can readily be extended 
to electron spin states in coupled quantum dots; 11 : 12 : 13 
interacting with a microwave cavity or an optical cavity. 

This work is organized as follows. In Sec. [U] we cal- 
culate the optical response of the coupled qubits in the 
dispersive regime, where the cavity mode is far off res- 
onance with all transitions taking place in the coupled- 
qubit system. We focus on the case of identical and of 
non-identical qubits with finite direct coupling. For iden- 
tical qubits we derive the fidelity for a cavity-mediated 
two-qubit gate in the presence of the direct qubit inter- 
action. In Sec. IIIII we discuss the control of the direct 
qubit interaction via the ac Stark shift for qubits with 
.9i 7^ <?2- In Sec. IIVI we discuss possible experimental 
implementations of Eq. (TlJ in physical systems and de- 
rive explicit expressions for the energy scales in Eq. (TlJ 
for quantum dots, quantum shells, and superconducting 
qubits in a cavity. We conclude in Sec. 
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II. SCHRIEFFER- WOLFF TRANSFORMATION 
OF THE HAMILTONIAN 

We consider the dispersive regime, where the cavity 
mode is non-resonant with all transitions of the qubit sys- 
tem and the coupling of cavity and qubits in Eq. (JTJ can 
be treated perturbatively. For a single qubit, the cavity 
resonance was shown£ to shift by ±<7 2 /A, where A is the 
photon detuning, depending on the state of the qubit. 
Measurement of the cavity resonance hence provides a 
readout mechanism for the qubit which has been demon- 
strated for Cooper pair boxes at the charge degeneracy 
points For two qubits which are not directly coupled to 
each other [J± = J z = in Eq. Q], but both of which 
are coupled to the cavity mode with g\ = g% = g ^ 0, 
the cavity-qubit coupling can be integrated out to lead- 
ing order by a Schrieffer- Wolff transformation U . The 
resulting effective Hamiltonian is£*& 



H' 
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with Ai t 2 = f2i ( 2 — u>. Equation implies that (i) 
photon emission of one qubit into the cavity and sub- 
sequent absorption by the other qubit gives rise to an ef- 
fective cavity-mediated inter-qubit coupling; and (ii) the 
cavity resonance is energetically shifted relative to u by 
^(-A^-A^ 1 ), ffW-Aa 1 ), ^(-A^+A^ 1 ), and 

lit), 
A 2 | > 



-A, 1 -A7 1 ), g 2 (A 1 1 ;>y 
g 2 (Aj" + Aj ) for the two-qubit states |tl), 



and | jf), respectively. In particular, for |Ai 
g 2 |(A 1 _1 + A^" 1 )!, spin flip-flop transitions between the 
qubits are weak and the cavity resonance can be used for 
the simultaneous readout of both qubits if the cavity loss 
rate is sufficiently small that all four frequencies can be 
resolved £ 



A. Identical qubits with direct coupling 

We now turn to the limit of identical qubits (Qi = 
Q2 = ^ and gi = 92 = <?) with finite direct coupling. 
Because 



N = ala + s lz + 



S2, 



(3) 



commutes with H, the problem can be solved exactly 
by numerical diagonalization of the Hamiltonian in the 
four-dimensional subspaces with given eigenvalue N of N 
(see below). First, we provide a perturbative analytical 
solution in the limit of weak coupling where the detuning 
of lu relative to all transitions of the coupled qubit system 
is large compared to gVN- For g = 0, H is diagonal 
in the sing let - triplet basis, \S) = (| U) - I it»/A 
\T-) = 1 ||), \T ) = (| tl) + I IT))/V2, and \T+) = | TT) 
with energies E s = -J*/4 - Ji/2, E- = -0 + J 2 /4, 
E = -J*/4 + Jl/2, and E + = Q + J z /A. We now 
generalize the Schrieffer- Wolff transformation of Refs.[fj||l 
for finite J± and J z by defining 



H' = UHU\ 



(4) 



where 



U = e J 



(5) 
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2A_ 



{a [s+(l - 8 2 ,,) + s+(l - «!,,)] - H.c.} 



+ 2X" {« + + * 2 + (l + *!,*)] - H.c.} 



with A = — lj, J — (Jj_ — J z )/2, and where A± = 
A =p J denotes the energy difference between the triplet 
state \T±) and \Tq). By construction, the Schrieffer- Wolff 
transformation removes the cavity coupling to leading 
order. Expanding H' to second order in g, we find 



H' 



Q(h + h, z + s 2 ,z) - An + (s+s 2 + s 1 §%) + J z h,zS2,z 



(6) 



n + 1 



-^A (Si,z + s 2 ,z) + 2g 2 {2n+l) (— - ] .s, . : .s. 



h + 1 
A_ 



n 

A+ 



'1 °2 ~ J l ^2 



1 1 

a7 _ a - 



(c^s+sj + a^ 2 s a s 2 ) 



where n — a^a is the photon number operator. The sec- 
ond line of Eq. © is diagonal in the singlet-triplet basis 
and can be interpreted as ac Stark shift and Lamb shift 
of the qubit states due to the cavity mode. The third 
line of Eq. © describes two-photon transitions between 



\T+) and |T_>. 

The expansion in g leading to Eq. © is valid as 
long as gVN < |A±|. Unless g 2 N\A^ - Al x | > 
|2A + ff 2 7V(A^ 1 + Al 1 )|, two-photon processes described 
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FIG. 1: (color online). Energy level spectrum in the subspace 
N — 1 as a function of A for g — 0.05Q and (a) J± = 0.2Q, 
J z = 0, and (b) J± = 0.3O, J z = 0, respectively. Solid lines 
indicate the exact energy eigenvalues, dashed lines the ap- 
proximate analytical values obtained with Eq. JSJ. The three 
anticrossings discussed in the text are indicated by arrows and 
circles in (b). 

by the last term in Eq. (jfJJl can be neglected. For a cavity 
mode detuning which is off-resonant with all bare tran- 
sition energies, Eq. © implies that the cavity resonance 
experiences a shift 



lu — > uj + 2g 2 x 






for |S), 


a; 1 


for \T + 




for |T ) 


-a: 1 


for |T_ 



Similarly to two qubits with different level spacing 
ill 7^ ^2 and without direct coupling, J± — J z = 0, 
see Eq. @, the cavity resonance splits into four lines, 
depending on the state of the coupled qubits. In stark 
contrast to qubits without a direct interaction, however, 
measurement of the cavity resonance projects the qubits 
onto the singlet-triplet basis rather than the product ba- 
sis of s z eigenstates. A resonance shift by 2g 2 (Al 1 — 
A^ 1 ), corresponding to the \T ) state, indicates that the 
coupled qubits are in a maximally entangled (and also 
super-radiant) state. In contrast, the singlet state \S) 
(being a sub-radiant state) decouples from the radiation 
field due to symmetry reasons and does not induce a shift 
of the cavity resonance line. However, such a proof of en- 
tanglement would be rather indirect because it relies only 
on the energy levels of the system rather than quantum 
state tomography. 

In order to verify the range of validity of our results, 
we compare the spectrum determined from Eq. (jSJ in the 
dispersive regime with exact diagonalization. In Fig. ^ 
we show the exact level spectrum (solid lines) as a func- 
tion of A for g = 0.05ft, J z = and (a) Jj_ = 0.2ft 



(a) (b) E[n] 




A [fl] 



FIG. 2: (color online), (a) Triplet level scheme in the "dressed 
state" picture, (b) Magnification of the exact anticrossing 
(solid lines) of |T_; 2) and |Tb; f) shown in Fig. 0(b), in com- 
parison with a two-level description (dashed lines) . 

and (b) J± = 0.3ft, respectively. The analytical results 
valid in the perturbative regime of large detuning from 
all resonances are indicated as dashed lines. 

As is evident from Fig. ^ the expansion in g breaks 
down if the cavity mode approaches resonance with one 
of the triplet transitions. On varying A, the spec- 
trum shows three anticrossings. Two strong anticross- 
ings with splitting g^j2(n + 1) (gy/2n) for A + J = 
(A — J = 0), respectively, corresponding to strong mix- 
ing of the states |T_;n + I) and \To;n) (\To;n) and 
|T+; n— 1)) by one-photon absorption. In addition, there 
is a weaker anticrossing with splitting 2g 2 y / n(n + 1)/J 
between |T_; n+1) and \T + ; n ~ I) at A = due to two- 
photon absorption processes described by the last term 
in Eq. © . In Fig. |2 we show all relevant transitions and 
the exact energy eigenvalues in comparison with approxi- 
mate expressions of a two-state model which are obtained 
from a diagonalization in a truncated basis, taking only 
the two nearly-degenerate states into account. 

In Fig. we illustrate the ac Stark shift as described by 
the second line of Eq. © as a function of the cavity pho- 
ton number n and the frequency lj. As two-photon tran- 
sitions are negligible, the energy eigenvalues only depend 
on the average photon number n. We note that for large 
photon numbers in Fig. |3 g^/n « A/2 is reached, for 
which the perturbative coupling to the cavity is not rigor- 
ously satisfied. We also show the ac Stark shift obtained 
for the optical interaction without the rotating-wave ap- 
proximation (RWA), which we derive in Appendix iBl 

B. Cavity-mediated two-qubit gate 

Before proceeding to the more general case of non- 
identical qubits coupled to a cavity, we discuss the effect 
of the direct qubit coupling on the proposed usage of the 
cavity as a quantum bus to couple the two qubits££ As 
shown in Ref. U a \J zSWAP two-qubit gate can be re- 
alized for J z = Jj_ =0, when the two qubits are tuned 
into resonance with each other (while maintaining off- 
resonant coupling with the cavity) for the interaction 
time i int = 7rA/4g 2 . However, in the presence of the 
direct qubit interaction, the two-qubit dynamics is de- 
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FIG. 3: Energies of two coupled identical qubits according to 
Eq. JSJ (solid lines) and for the optical interaction without 
the RWA (dashed lines), see Eq. 1B7L of the triplet states 
(a) \T-), (b) \T+), and (c,d) \To) as a function of the photon 
number n or the frequency ui of the cavity mode. The dotted 
lines indicate the energy levels in the absence of the coupling 
to the cavity, respectively. The parameters used are J = 
0.001ft and g = 0.0002a Further, uj = 0.96ft in (a,b,c), and 
n = 10 4 in (d). 



termined by Eq. @) rather than Eq. J2J), and the effec- 
tive interaction via the cavity is accompanied by the di- 
rect qubit coupling. For simplicity, we focus on the case 
J z = here, which is realized, e.g., for capacitively cou- 
pled Cooper pair boxes, see Sec. lIVCl The total effective 
two-qubit coupling energy in Eq. 10 is given by 



h(=(n Iff'UT) 



A 



hJ ± 



l 



A_ 



n 

A~ 



(8) 



l-^(2n + l) 



where the second line is valid for J± <C 2A. We quan- 
tify the fidelity of the ViSWAP gate proposed in Ref. @ 
as /(i) = <TI |^_C/(t)| U), where 



iSWAP 



is the 



ViSWAP 

ideal gate operator and U{t) is the time-evolution opera- 
tor for the two qubits due to Eq. ©. Using the first line 
of Eq. (JSJ, we obtain for the gate fidelity 



/(tint) 



COS (C^int 

r 2 7 - 



tt/4) 



(9) 



J 2 A 2 



128g 4 



1 



1) 



Here, we also give the asymptotics for J± <C 2g 2 / A, when 
the cavity-mediated interqubit coupling dominates over 
the direct exchange-like coupling quantified by J±. We 
notice that /(tint) ~ for Jj_ w 16g 2 /A. For exam- 
ple, for a coupling HJ± = 0.1 neV (hJ± = 0.2 neV) and 
the parameters^ hg ss 0.01 /ieV, HA w 1 //eV, we ob- 
tain /(t int ) ~ 0.924 [/(i int ) « 0.708] for a photon num- 
ber n = 10, indicating that a residual direct inter-qubit 



coupling may lead to gate errors in cavity-mediated two- 
qubit gates. However, this error may be significantly re- 
duced by adjusting ij nt , as discussed in the following. 

In order to reduce the fidelity loss discussed above, a 
two-qubit gate could be implemented based on the to- 
tal two-qubit coupling energy Eq. JSJl, in particular for 
systems for which the qubit interaction via the direct 
coupling and via the cavity are of comparable strength. 
If both qubit energies are strongly detuned from ui and, 
additionally, | $~2x — O2 1 3> J±, the direct qubit interaction 
is strongly suppressed along with the effective interaction 
via the cavity. Starting in such a situation, a v 7 iSWAP 
two-qubit gate could be realized with Eq. © by tuning 
the two qubits into resonance and establishing simultane- 
ously a smaller detuning of the qubits to the cavity. For 
the total two-qubit interaction, we readily obtain the re- 
quired interaction time for a V iSWAP gate (up to phase 
factors) along the lines of Ref. y as t[ nt = 7r/4£. If a non- 
zero coupling J z is introduced, the block-diagonal form 
of U{t) is preserved, and the 2x2 subspace spanned 
by I Tl) an d I IT) contains a time-dependent phase factor 
exp (— iJ z t/2) with respect to the 2x2 subspace spanned 
by I ft) and | J.J.). We note that for the special case 
Jz/4: = J±/2 + g 2 [(n + 1)/A_ - rc/A+], with a fixed 
photon number n, the total coupling of the two qubits 
is of the Heisenberg form, enabling a VSWAP two-qubit 
operation rather than -\ASWAP. 



C. Non-identical qubits with direct coupling 

Here, we provide the Schrieffer- Wolff transformation 
for the general case of arbitrary system parameters f2i, 
gii and <?2 in the Hamiltonian Eq. For vanishing 
qubit-cavity coupling, gi = g% = 0, the eigenstates of the 
coupled qubit system are 



\T+) = ITT), 

\T-) = III), 

= cos(0/2)|T|)+sin(0/2)||T), and 

\1> a ) = sin(0/2)|T|)-cos(0/2)||T) 



(10) 



with eigenenergies 



hE ± = h^±h^±^. 
4 2 



(11) 



h,E s/a = -h^±n^(n 1 -n 2 ) 



respectively, where tan</> = Jj_/(f2i — Q2)- 

In contrast to the symmetric case discussed in Sec. [HJ 
where the singlet state decouples completely from the 
cavity field and only two transitions are allowed by se- 
lection rules, there are four allowed transitions |TL) 
\ips/a) an d |Vty<») l^+) in the general case, see Fig. 0J If 
the cavity mode is off resonance with all four transitions, 



(£4 



E. 



/a) 



0J- 



(E s/a -E-)\ > giVn,g 2 y 



with n the characteristic photon number in the cavity, 
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the cavity mode can be treated perturbatively. In the 
eigenbasis Eq. (|1(J|) . the Hamiltonian Eq. reads 

H = J2 hE i\ T i)( T i\+ J2 KEi^&il+tiLjtfa (12) 

i— ± i— s,a 

+ah ( 7l |^ s ) (T_ | + 73IVO (T_ | + 7 3 |T+) (V> s | 
+74|r + )(^ a |)+a t ft(7i|r_)(V s |+72|T-)(^ a | 

+73|^)(T + |+ 7 4|Va)(T+|), 



E 



with the coupling constants 



7i.4 = ±cos(^/2)g! +sin(0/2)g 2) 
72,3 = sin(0/2)5! =fcos(0/2)5 2 . 



(13) 
(14) 



To first order in 7 , , the Hamiltonian is diagonalized by 
the Schrieffer- Wolff transformation 



H' 



aB-a* B f jj -aB+o t 



with 



7l |^)(T_|+7 2 |^)(T_| 
+7 3 |T + )(^|+74|T + )(^ a | 

The coefficients 7 j are defined by 



(15) 
(16) 
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72 
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(E a 


-E-) 


" A a _ 


73 = 




73 




_ 73 


UJ — 


(E + 


-E s ) 


~ A +s 


74 = 




74 




_ 74 


UJ — 


(E + 


- E a ) 


~ A +a 



(17) 



By assumption, |7j| <C 1 because of the off-resonance 
condition. The expression for H' can then be truncated, 
neglecting terms of third and higher order in 7$. Up to 
second order in the coupling coefficients 7 j, we obtain 



H — H + i?s,L + E[ 2 - P h + H n 



(18) 



Here. 



H = J2 ^E % \T % )(T % \ + J2 m\il>i){ik\ + (19) 



is the unperturbed Hamiltonian for vanishing cavity cou- 
pling. 

#S,L = S (o)a + 1) (7373 + 7474) \T+){T+\ 

-M t a( 7 i 7 i+7272)|T_){T_| (20) 
+h [a f a (7171 - 7373) + 71 71] \ip a )(ipa I 
+h [a f a (7272 - 7474) + 7272] \ipa}(ip a \ 
denotes the Stark and Lamb shift of the four basis states. 



2-ph 



fi 737i ~ 7i73 + 7472 ~ 7274 / 21 s 
x(a 2 |T + )(T_|+at2| T _ )(r+ |) 
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FIG. 4: Energy level scheme of two coupled qubits for the 
general case. The arrows indicate transitions under photon 
absorption and emission. Contrary to the symmetric case 
discussed in Sec.[n] for the general case there are four allowed 
transitions with coupling constants 7$, i = 1 . . .4. 



describes two-photon transitions between \T + ) and |T_) 
(a^a + 1) (7172 + 7271) - a f a (7374 + 7473) 



X(|VaW«l + IV'a><^|) 



(22) 



is a term which gives rise to transitions between the ba- 
sis states \if> 3 ) and \ip a )- Its physical meaning will be 
explained in more detail in Appendix ^ We note that 
for J± = the effective interaction of | f J.) and | if) via 
the cavity is given by H m { x . On a technical level, -ffs,L 
and i?mi x correspond to virtual processes corresponding 
to the emission and re-absorption of a photon. 

The cavity line shift induced by Hs,l for the coupled- 
qubit eigenstates is given by 

( 7f/A +s + 7 J/A +a 
™ 7l 2 /A s _ - 7 f/A a _ 
7l 2 /A s _ - 7 3 2 /A +s 
7 2 /A Q _ - 7 4 2 /A+a 



UJ + < 



for \T+) 

for |T_) 

for \ip s ) 

for \ip a ) 



(23) 



where the coupling coefficients 7$ are given by Eqs. i|13|) 
and 1)14(1 . and the effective detunings A nm are de- 
fined by Eq. ((17(1 . However, above states are only 
eigenstates of the qubits coupled to the cavity when 
H m i x + H 2 - P h can be neglected. For \J Z \ < \J n \ = 

(1/2)1^/(01 - n 2 ) 2 + Ji 1 < |A| = |(ni+n a )/2— w|, the 

mixing energy of |^ s ) and \ip a ) due to i? m i x is approxi- 
mately given by h[sm(j)(gf — g 2 )(2fi+l) — 2gig 2 cos</>]/2A. 
It is typically small due to the perturbative coupling con- 
dition and vanishes for coupled identical qubits (</> = 7r/2 
and gi — g 2 ). Similarly, in the limit | JoJ <C |A|, the mix- 
ing energy of |T + ) and |T_) due to H 2 _ p h is proportional 
to hsin(j)(gf + g 2 )Jn/A 2 and is therefore negligible. 

The above analysis has been focussed onto coherent 
quantum dynamics. In real systems, however, the rate T 
for spontaneous photon emission and the decay rate k of 
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the cavity field lead to spectral line broadening. For the 
Jaynes-Cummings model in the dispersive regime, the 
decay rates of the two dressed states of lowest energy 
(i.e., in the subspace of one energy quantum) are already 
given in Ref. |5|. Similar to the readout of two qubits 
with no direct interaction^ we expect that the line shifts 
for directly coupled qubits can be well resolved if they 
are larger than k + T. A more detailed treatment of 
dissipation is beyond the scope of the present work. 



III. CONTROL OF THE DIRECT QUBIT 
INTERACTION VIA PHOTONS 

As a second application of two qubits interacting with 
a cavity, in addition to the readout discussed above, we 
show that the qubits can be tuned into and out of res- 
onance by controlling the photon number in the cav- 
ity mode if the qubit-cavity couplings are not identical. 
Here, we focus on the case of identical qubit level spac- 
ing, rii = f^2 with strongly asymmetric cavity couplings, 
gi = g ^ 0, gi = 0. As shown in Appendix lAl the Stark 
and Lamb shift and mixing Hamiltonian, fls,L + ^mix, 
then induce a relative energy shift of the states | f|) and 
| J,f) which can be used to control the level admixing 
induced by the direct coupling J±. 

For illustration, we focus on a cavity mode close to 
resonance with the transitions from \T + ) to | Jt), \u> — 

+ J z /2)\ < \uj- (fi- J z /2)\, but still in the dispersive 
regime, \g^/n/[uj — (Q + J z /2)] \ <§; 1. Then, the ac Stark 
shift due to qubit 1 is given by 



9 



to - (0 + Jz/2) 
x [ata||T)UT |-(6 f fi + l)|T + )<T + |]. (24) 



If the number of photons in the cavity is n such that 

-,2 



9 



uj-(n + j z /2) 



> J 



(25) 



the two qubits are effectively decoupled due to the ac 
Stark shift, since level mixing due to J±_ is strongly 
suppressed. Extracting the cavity photons on a time 
scale short compared to r sw = 2ir/J± and maintaining 
the empty cavity state for a predetermined time ensures 
that, in the absence of an ac Stark shift, the direct in- 
tcrqubit coupling J± induces quantum coherent oscilla- 
tions between | fJ.) and | J,|). By reinserting a number 
of photons satisfying Eq. (|25J) faster than the time r sw 
into the cavity, the two qubits can be decoupled again. 
Based on the latter inequality we define a minimum pho- 
ton number for optical switching of the direct coupling, 
"min = J±\u} — + J z /2)\/g 2 . However, in order for the 
analysis leading to Eq. (|2~4f> to remain valid, the perturba- 
tive coupling condition n m i n g 2 <C \u> — (il + J z /2)\ 2 must 
be satisfied. From this follows |cj — (fi + J z /2)| ^> J± 
and, using the definition of n m i n , n m i n ^ (J±/g) 2 . Ob- 
viously, decoupling of the two qubits requires a larger 





QD 


QS 


CPB 


T sw [ps] 


0.4. . .4 


4 ... 40 


400 .. . 4000 




10 3 . . . 10 4 


10 . . . 10 2 


10 5 . . . 10 6 


hJ ± [meV] 


1 ... 10 


0.1... 1 


10~ 3 . . . 10~ 2 


hg [meV] 


0.1 


0.1 


lO" 5 


%-(fi + J z /2)\ [meV] 


10 


1 


0.01 



TABLE I: Switching time r sw = 2ir/ J± and lower bound (see 
text) for the minimum photon number J±\lu - (fi + 

J z /2)\/g 2 to decouple two qubits with fii = Q.2 and <?i = g 7^ 
(/2 = via the ac Stark shift in a cavity. The considered types 
of qubits are excitons in coupled quantum dots (QD), excitons 
in quantum shells (QS), and capacitively coupled Cooper pair 
boxes (CPB) with estimates for J± and J z from Sec. IIVI The 
estimates for kg are based on Refs. 1151181 



photon number n > n m i n in the case when both coupling 
constants are non-zero (but still satisfy \gi\ ^ |<?2|)- I n 
Section ITVI we study three example systems for the cou- 
pled qubits: excitons in quantum dots, excitons in quan- 
tum shells, and capacitively coupled Cooper-pair boxes. 
For those systems, we calculate the parameters J± and 
J z for typical coupling mechanisms. Based on these esti- 
mates, we present in Table [I] estimates for the switching 
time r sw and the minimum photon number n m j n required 
to switch the direct coupling for the energy scales of h J± 
obtained in Section llVl In Tabled we consider detunings 
satisfying \u> — (Q + J z /2)\ > J± within the given range 
of Ji , respectively. For the maximum values of J± , these 
detunings and rt m in do not strictly satisfy the condition 
for perturbative coupling to the cavity and thus provide 
a lower bound for n m j n . 

While the addition of photons to the cavity mode can 
be realized via external irradiation, a fast extraction re- 
quires a cavity design enabling a quick reduction of the 
Q-factor. For optical microcavities containing Bragg mir- 
rors, this could be achieved, e.g., via optical switching 
of the reflectivity^ Such a scheme would facilitate an 
all-optical control of the interaction of two qubits with 
possible applications as an optically controlled two-qubit 
gate. However, a high degree of control on the parame- 
ters of a high-Q cavity is required to perform such a gate 
operation via a cavity mode. In particular, the quick 
reduction of the Q-factor still needs to be demonstrated 
experimentally. As an interesting alternative, an ac Stark 
shift for exciton qubits in quantum dots or shells can also 
be induced by a laser instead of a cavity. The switching 
"on" and "off" of the direct qubit coupling might then 
be accomplished with ultrafast optical techniques. As a 
recent first step into this direction, the laser-induced ac 
Stark effect has been demonstrated experimentally for 
single GaAs quantum dotsiii 

Complementary to the decoupling of two identical 
qubits with asymmetric coupling to a cavity (or laser) 
mode, the ac Stark shift could be exploited to bring two 
non-identical qubits into resonance, as already suggested 
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by Nazir et ali^ for excitons in quantum dots coupled 
to a laser. The control of the direct interaction J± could 
then be achieved via the shaping of microwave or laser 
pulses, for superconducting or exciton qubits, respec- 
tively. We note that in our work the exact Schrieffer- 
Wolff transformation is provided for the coupled system 
with and without the RWA for the optical interaction, 
while Nazir et al. provided a numerical solution within 
the RWA. 

The coupling constants g\ and gi of the qubits to the 
cavity may be intrinsically different due to fluctuations in 
the composition, size, and geometry of the structure con- 
taining the qubit, which may randomly occur or might 
be induced on purpose in the fabrication process. In 
principle, asymmetric coupling of the two qubits to the 
cavity could also be achieved due to the different qubit 
locations with respect to the cavity field. Still, the dis- 
tance between coupled quantum dots or superconducting 
islands is typically much smaller than the cavity-photon 
wavelength, which complicates this approach. A strongly 
asymmetric cavity coupling could also be achieved via 
the misalignment of the electric dipole moment of one 
qubit relative to the electric field in the cavity, e.g., for 
Cooper-pair boxes in a line resonator. Ideally, (say) d2 
is oriented perpendicular and di is oriented parallel to 
the cavity field, such that the coupling of electric dipole 
transitions to the cavity is strongly asymmetric. 



IV. MICROSCOPIC MODELS 

We next derive explicit expressions for the characteris- 
tic energy scales in Eq. Q for qubits defined in terms of 
exciton states of a given spin in quantum dots or quantum 
shells and for superconducting qubits at the charge de- 
generacy point. While most quantum computing schemes 
for quantum dots rely on the electron spin as qubit j2S we 
focus on charge-based qubits defined by the presence or 
absence of an exciton (see, e.g., Ref. [21J), because such 
states couple to the cavity mode by photon emission and 
absorption. Here, we rely on the fact that the exciton 
spin relaxation time is typically much longer than the 
exciton lifetime in quantum dots at low temperatures 
and low magnetic fields^ Optical selection rules then 
ensure the coupling of the polarized excitons to a cav- 
ity mode with a suitable polarization, see, e.g., Ref. l23l 
for typical optical transition matrix elements for quan- 
tum dots. The integration of quantum dots into high- 
Q cavities is still a considerable challenge, but has been 
demonstrated for self-assembled quantum dots^ and also 
appears to be feasible by deposition of chemically synthe- 
sized nanocrystals onto defect modes of photonic crys- 
tals^ For definiteness, we focus on the latter system in 
the following. 

The analysis can be easily extended to electron spin 
qubits in quantum dots, for which the electron exchange 
interaction leads to an effective coupling of the Heisen- 
berg form [i.e., J± = J z /2 in Eq. JJJ], provided that the 



single-dot charging energy is larger than the tunnel ma- 
trix element, which has already been studied in detail in 
Refs. ITTIIT^fia . For electron spins, a static magnetic field 
gives rise to Zeeman splittings Mli , typically on the order 
of tens of fj,eV for fields on the order of 1 Tesla, and the 
two spins couple, e.g., via magnetic dipole transitions to a 
microwave cavity. Alternatively, cavity-mediated spin in- 
teractions are implemented by two-photon Raman tran- 
sitions via trion states^ We note that the spatial separa- 
tion of coupled dots is typically too small to enable the 
application of individual laser fields to the two dots. Still, 
effective spin interactions could be mediated via the cav- 
ity if the electrons in the two different dots have different 
Zeeman splittings, e.g., due to magnetic dopants^ (giving 
different gyromagnetic factors), magnetic field gradients, 
or inhomogeneous hyperfine polarizations, since then two 
laser fields with different frequencies could be applied to 
the two dots simultaneously to establish the Raman tran- 
sition conditions. For the special case of identical dots, 
these conditions can be satisfied for both dots even by 
the same laser field. 



A. Quantum dots 

The confinement of electrons in nanocrystal quantum 
dots gives rise to well-defined quantum size levels with 
a spacing which is typically of order 0.1 eV for the con- 
duction band. We first consider pairs of nanocrystals 
interacting with a o~ circularly polarized cavity mode 
[Fig. OJa)], assuming that the symmetry axis of both 
nanocrystals is aligned with the propagation direction of 
the cavity photon. For definiteness, we will focus on ma- 
terials with valence band states described by the spherical 
Luttingcr Hamiltonian, which is a good approximation 
for, e.g., CdSe. The cavity mode is assumed to be suffi- 
ciently close to the absorption edge that only the lowest 
exciton transition from the 1S3/2 valence band multiplet 
to the conduction band ground state lS e must be taken 
into account. Crystal or shape anisotropy are assumed 
to split the light-hole (lh) and heavy-hole (hh) states at 
the r point and only transitions involving hh states are 
considered. The ground state of the nanocrystal, i.e., the 
exciton vacuum, is denoted as 1 1), while the one-exciton 
state is denoted as | \). 

The transverse coupling term in Eq. corresponds 
to coherent exciton transfer between the nanocrystals. 
Possible microscopic mechanisms are direct transfer of 
electron and hole, e.g., via a coupling molecule^*2Si2i 
or resonant Forster transfer In order for Eq. 

to capture the essential features of molecularly coupled 
nanocrystals, indirect exciton states must be negligible. 
For pairs of identical nanocrystals this implies that the 
transfer matrix elements for both the conduction and va- 
lence band, t c and t v , respectively, must be small com- 
pared to the exciton binding energy U ~ 1.6e 2 /47reoea, 
with the nanocrystal radius a. Data for CdSe nanocrys- 
tals coupled by benzene rings2£ are consistent with t c ~ 
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FIG. 5: (color online). Schemes for coupled qubits: (a) Pair 
of quantum dots interacting with one cavity mode as dis- 
cussed in Sec. II V Al (b) Nested quantum shells as discussed 
in Sec. IIV Bl (c) Cross-section through nanocrystal contain- 
ing a pair of nested quantum shells, (d) Two Cooper-pair 
boxes in a line resonator as discussed in Sec. IIV Cl 

80 meV while U ~ 160 meV for a = 1 ran. Because t c 
and U are of the same order of magnitude, the transfer 
matrix element for an exciton, HJ± = 2t c t v /U ~ t v , is 
determined by t v . While t v is expected to be smaller 
than t c because of the larger hole mass, atomistic calcu- 
lations for molecularly coupled nanocrystals have shown 
that electron and hole transmission probabilities are 
comparable^ The exciton transfer energy HJ± can thus 
presumably reach values of order 10 meV. 

Alternative mechanisms for exciton transfer are 
Forster processes which result from the electrostatic 
dipole interaction HamiltonianSSiSLSi 

( ri )(r 2 )--l(( ri ).R)(( r2 ).R) 



hJ± = 



AneaeR 3 



(26) 

where R denotes the relative position of the nanocrystals 
and (ri ^) are the dipole transition matrix elements be- 
tween the exciton vacuum and the one-exciton states in 
nanocrystals 1 and 2, respectively. Because the symme- 
try axes of both nanocrystals are aligned by assumption 
[Fig. [Sf a)], the Forster transfer in Eq. (|26|l conserves the 
exciton spin^Si The dipole matrix element connecting 
the —3/2 state of the IS3/2 valence band multiplet with 
the s z = —1/2 conduction band ground state is 

( 1 

(15 e ;T|r|15 3/2 ;-3/2) = (S\x\X)- ' 



1 

71 



In 



'3/2 



V 



is, 



(27) 



denotes the interband matrix element. The or- 
bital overlap of the envelope wave functions of lS e and 
the L = component of the IS3/2 multiplet, 



S , 3/2;iS' e 



= / drr 2 ip 1Sc {r)Ro{r), 



(28) 



is typically of order unity. For CdSe nanocrystals with a 
high-bandgap capping layer, I\s 3/2 -is c = 0.88i£i Insert- 
ma; Eq. (J23 into Eq. ®, we find for the Forster transfer 
energy 



e 2 \(S\x\X)f 
4:Tre eR 3 



(29) 



For CdSe, |(iS'|a;|X)| ~ lnm, such that hJ± = 
1.44eV/e(i?[nm]) 3 is of order lmeV for R = 5nm and 

6 = 9. 

The longitudinal coupling term oc HJ Z in Eq. Q repre- 
sents the biexciton shift. Dominant contributions are the 
electrostatic quadrupolar interaction and van-der-Waals 
interactions which correspond to simultaneous virtual 
transitions between the IS3/2 and IP3/2 valence band 
multiplets in both dots. The electrostatic interaction re- 
sults from the charge density of the 153/2 valence band 
state,2ii22 

\1S 3/2 :F Z =3/2) (30) 
= R (r)\L = 0, j = 3/2, F = 3/2, F z = 3/2) 
+R 2 {r)\L =2,j = 3/2, F = 3/2, F z = 3/2). 

i?o,2 are the radial wave functions for an orbital angular 
momentum L = 0, 2 of the envelope wave functioni^SiSi 
The electrostatic quadrupolc moment for a spherical 
quantum dot with radius a and a high-bandgap cap is 



drr 4 R (r)R 2 (r) 



(31) 



for vanishing boundary conditions. For a spherical CdSe 
nanocrystal with radius a, Q zz — — 0.06ea 2 . The biexci- 
ton shift resulting from quadrupolar electrostatic inter- 
actions^ is for a <C R given by 



hJ z 



(Q«/e[mn 2 ]) 2 
(i?[nm]) s 



8ire Q R 5 



2.16 eV 



(32) 



Because the radial integral in Eq. I|31|) is typically small, 
the electrostatic contribution to the biexciton shift is 
smaller than 1 meV for characteristic parameters. For 
CdSe nanocrystals, a = lnm (2nm), and R = 5nm, 
HJ Z = 2.4 LieV (39.8 ^eV). However, Eq. is no longer 
rigorously applicable in this parameter range and may 
only provide an estimate for the characteristic energy 
scale. 

For these or comparable parameters, the dominant 
contribution to the biexciton shift is caused by inter-dot 
van-der-Waals interactions. Because the level spacing for 
valence band states is small compared to the conduction 
band, the dominant contribution involves simultaneous 
transitions between the IS3/2 and IP3/2 valence band 
states in both quantum dots. In second order perturba- 
tion theory in the hole-hole Coulomb interaction, 
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hJ 



z,vW 



. 47re ei? 3 , 



E 



\ 1 {lP 3/2 -F z M lP 3/2\F z ,2\x 1 X2 + yiy2-2z 1 z 2 \lS z/2 -i/2) 1 \lS V2 -?,/2) 2 \ 

AEi + AE 2 



(33) 



The matrix element is finite only for F Zt \ = F z>2 = 3/2, 
i.e., for processes which describe simultaneous transitions 
from |lS , 3/ 2 ; 3/2) to |LP 3/2 ;3/2) in both dots. AEi. 2 de- 
notes the energy difference between the corresponding 
hole levels in quantum dot 1 and 2, respectively. For a 
pair of identical dots, 



W 



2|<lP 3/2 ;3/2|z|lSg /2 ;3/2>| 



The transition matrix element, 



AE 



(34) 



lP 3 /a;3/2|*| 



lS 3/2 ;3/2) ~ £ drr 3 Ro(r)Ri 



(r) 
(35) 

depends on the specific material parameters, but is typi- 
cally of order 0.1a. For CdSe, (lP 3/2 ; 3/2|z|lS 3/2 ; 3/2) = 
0.17a, while AE = -27 mcV/(a[nm]) 2 . For a = 1 nm 
(2 nm) and R = 5 nm, the van-der-Waals contribution 
to the biexciton shift is hJ z , v w — —8.2fieV (— 132/xeV), 
larger than the electrostatic contribution derived above. 

We conclude that for nanocrystals with radius of order 
1 — 2 nm and a distance of order 5 nm, the energy scale 
for Forster transfer, 1 — lOmeV, is large compared to 
the inter-dot biexciton shift which is typically of order 
- 10 fxeV. 



B. Nested quantum shells 

Similarly to quantum dots, quantum-dot quantum 
wells or 'quantum shells' 33,35,36 provide quantum confine- 
ment for electrons and holes which allows one to define 
a qubit by the presence or absence of an exciton with 
a given polarization. While experimental studies have 
so far focussed on individual shells, with a low-bandgap 
material sandwiched between a high-bandgap core and 
cap, synthesis of nested shells is possible. Such a sys- 
tem is expected to exhibit properties similar to a pair 
of quantum dots [Figs. EJb),(c)]. Here, we derive the 
characteristic energy scales for exciton transfer between 
shells (HJ±) and the biexciton shift (hJ z ). All energy 
scales typically depend on the inner and outer radius of 
a shell, the potential barrier, and the band masses inside 
the shell and the barrier. For definiteness, we consider 
shells with width Wi small compared to their inner radius 
a,, where i = 1,2, in the limit of high potential barriers 
where both conduction and valence band states are well 
confined inside the shells. Then, the wave functions of 



the lowest conduction band state is 

, , s nrr, — r sin[7r(r — aA/wA T , ... 
^•(r) = V(2M) ^^Foo(r), 



(36) 



while the wave function of the 15' 3 / 2 valence band states 
with total angular momentum quantum numbers F = 
3/2 and F z is 



a IT\ 77. \ sin i n ( r - a i)/ w i 



(37) 



x (|L = 0;j = 3/2;F = 3/2;F z ) 
+\L = 2;j = 3/2;F = 3/2;F z }). 

The energy splitting to the lP 3 / 2 multiplet is ft 2 (71 — 
2j)/moaf, where mo is the free electron mass and 71 and 
7 denote the Luttinger parameters of the shell material. 

Dominant mechanisms of exciton transfer between the 
shells arc electron and hole tunneling or Forster trans- 
fer. For a finite potential barrier Vo,i of width w, the 
tunneling matrix elements for electrons (i = c) and 
holes (i = v) are exponentially suppressed with a factor 
exp (— 2wy/2rrii(Vo,i — Eij/H). For realistic structures, 
such as CdSe/ZnS/CdSe nanocrystal heterostructures, 
the electron tunneling matrix element can be of order 
~ 10 — 50 meVM' In order to calculate the matrix ele- 
ment for the Forster process, we invoke the identity 



1 



r x -r 2 



E 

l:rn 



47T 
2T+1 



^-m(ri)F /m (r 2 )- 



,z+i ' 



(38) 



where T\ < r 2 , i.e., the index 1 denotes the inner shell. 
Similarly to quantum dots, Forster transfer for nested 
quantum shells conserves the exciton spin. In order to 
estimate the characteristic energy scale of the Forster 
matrix element, we calculate the transition matrix ele- 
ment for the dominant 1 = 1 — term in Eq. (|38|l and ap- 
proximate the hole state by its L — envelope function 
component. Then, 

1 4-TT 1 2 

(iuta I- — z-rl Ti;| 2 ) - ^r- \(S\x\x)\* 



|ri - r 2 \ 

/•a 1 +wi 



3 3n w\w 2 



(39) 



X / drsin [7r(r — ai)/ti>i] 

x r +M ynr 2 Kr-a 2 )M] 
/ r 3 

^2\[S\x\X)l 
9a 3 

The approximation in the last line in Eq. I|39|l is valid in 
the limit of narrow shells considered here. (5|x|X) de- 
notes the Kane interband matrix element. Assuming that 
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the dielectric constant of shells and potential barriers are 
comparable, Eq. I|39[) implies 



2e 2 \(S\x\X)f 
36-7T eoe a| 



(40) 



which only depends on the radius of the outer shell 
but is independent of the size of the inner shell. For 
CdSe, e = 9 and |(5|x|X)| ~ lnm, such that hJ±_ ~ 
35 meV/(a2[nm]) 3 . For small nanocrystals, where, e.g., 
0,2 — 4 nm, the Forster matrix element is of order h J± ~ 
0.5 meV. 

Similarly to quantum dots, the dominant contribu- 
tions to the biexciton shift for quantum shells are the 
electrostatic quadrupolar interaction caused by the non- 
spherical charge density of the holes and van-der-Waals 
interactions. With the expansion Eq. 1)38)1 ■ after integra- 
tion over the angular variables, we obtain for the electro- 
static quadrupolar hole interaction 



HJ Z 



47ree 25 J ai 



Ol+U>l 



a 2 +w 2 



dr< 



dr 1 r 4 1 R {ri)R2(ri) (41) 
Rp{r 2 )R 2 {r 2 ) 



47ree 25a| 



where, in the last line, the radial integrals have been 
evaluated for the limit of narrow shells considered here, 
Eq. (|37)l . and it is understood that the hole functions 
i?o,2 m Eq. (|41|l denote the corresponding components of 
the hole wave functions in shells 1 and 2, respectively. 
For ai = 2nm, a 2 = 3nm, and e = 9, ftJ 2 ,Q = 0.9 meV. 

The van-der-Waals interaction caused by simultaneous 
transitions between IS3/2 and IP3/2 hole multiplets in 
both shells can be calculated similarly to Eq. (|33|) in sec- 
ond order perturbation theory in the hole-hole Coulomb 
interaction. In terms of the radial wave functions for the 
IS3/2 and IP3/2 valence band multiplets, 



(— 

■ai+wi 



1 



/ dn r\ 



(L2+W2 

dr 2 

2 



3 / RoR 



AEi + AE 2 
1 2R 2 Ri 



2 5 
RqRi 2R 2 Ri 



3R 2 Rs 
10 



(42) 
2 



?>R 2 Ri 
10 



n 2 



36ai 



1 



~ ("47r<„< Vi:>«-, ) AE, - AE 2 ' 

The last line is valid in the limit of narrow shells, where, 
again, it is understood that the hole wave functions de- 
note the wave functions for the respective shells. For 
a% = 2 nm, a 2 = 3 nm, the characteristic hole level spac- 
ing is of order lOmeV. For e = 9, the characteristic 
energy scale for \hJ z ^ vW \ is lOmeV, significantly larger 
than the electrostatic quadrupolar energy derived above. 



C. Superconducting charge qubits 

Extending the model of Ref. 0, we assume that two 
small superconducting islands in the charge regime (so- 
called Cooper-pair boxes) are located at the center of a 
resonator and are separated by a distance r, see Fig. 
(d). In the charge basis, the total capacitive coupling of 
the two singly charged Cooper-pair boxes is given by a 
z-z coupling term with energyS*^ 



(2e) 2 C c 



c c )(4 2) 



c c )-ci 



(43) 



Here, C c is the capacitance that couples the two islands 



and Cj^ = Cji +Cj l> is the total capacitance of the single 
superconducting island i = 1,2, where C^f 1 (C^) is the 
capacitance of the island and the corresponding Joseph- 
son junction (the resonator) From the recent experi- 
ments using a single Cooper-pair box in a resonatoriii^ 
we infer ~ e 2 /40/ieV ~ 4 fF. We give a rough esti- 
mate for C'c in the following and approximate the Cooper 
pair boxes as two rectangular plates of length I and width 
w. The electric field along the perpendicular bisector of 
a wire of length I with charge density a is calculated as 
Ewirc(^) = ul/AneeQX^ x 2 + ^ 2 /4, where x is the distance 
from the wire. Integrating over the width w we obtain 



Epl a te(x) 



Q 



2tte£qIw 
— arcsinh 



arcsinh [ ± 
2x 



I 



2(x + w) 



(44) 



for the central electric field at a distance x from the bor- 
der of the plate [see Fig. [S] (d)], where Q = alw is the 
charge on the island. We obtain for the capacity of the 
two islands 



2 I dxE p i atc {x) 

2lT£EoW 



(45) 



where we use K(x) — x arcsinh(l/x) + arcsinh(x). For 
I = 0.1 /im, w = 1 /im, and r = 0.5 /im, we thus obtain 
C c w 15. 5e aF, where e should be understood as the 
average dielectric constant defined for the space between 
the Cooper pair boxes. For e = 9, C c w 0.14 fF. 

Using the qubit eigenstates at the charge degeneracy 
point, |t), ||), the capacitive interqubit coupling trans- 
lates into hJ x s\ x s 2 x with h,J x = E cc [Eq. This 
interaction mixes the states | tt) an d | II) with the 
mixing angle arctan( J x /[2(f2i + ^2)])- Since the single- 
qubit splittings are given by the Joscphson energy 
(» 30 ^teV in Ref.0), J x < 2(fii + fi 2 ) and the interqubit 
coupling maybe approximated by (HJ x /4:)(s^ §2 s^)- 
Therefore, capacitively coupled qubits interacting with a 
cavity may be mapped onto Eq. |{T{|. wherein J± — J x /2 
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and J z = 0. For C c « 0.14 fF and characteristic parame- 
ters for Cooper pair boxes in a resonator we obtain for 
the coupling energy HJ± — E cc /2 w 2.6 /j,eV. 



V. CONCLUSION 

We have studied two directly coupled qubits which 
interact perturbatively with a single off-resonant cavity 
mode. Based on the Hamiltonian Eq. JTJ we have pro- 
vided a diagonalization by means of a Schrieffer- Wolff 
transformation for various cases in order to obtain an ef- 
fective Hamiltonian describing the perturbative interac- 
tion of the cavity with the two-qubit system. For identi- 
cal qubits, the effective Hamiltonian is given by Eq. JJjJ, 
while for the more general case of non-identical qubits 
the effective Hamiltonian is given by Eqs. (|18f) - (|22[1 . 
Beyond the RWA for the optical interaction, the effec- 
tive Hamiltonian for non-identical qubits has been ob- 
tained as Eqs. I|B6(I - l|B9(l . Based on these results, we 
have shown that the state of the coupled qubits can be 
read out from the photon dispersion, the coupled qubit 
state being projected onto a singlet-triplet basis. Fur- 
ther, we have investigated the eff ect of th e direct qubit 
interaction on a cavity-mediated v zSWAP gate and have 
shown that a two-qubit gate can also be realized with 
the total qubit-qubit interaction for a finite residual di- 
rect inter-qubit coupling which might be present for, e.g., 
Cooper-pair boxes in a resonator. For asymmetric cou- 
pling of the qubits to the cavity, we have discussed the 
control of the direct qubit interaction via the ac Stark 
effect. Finally, we have derived explicit expressions for 
the energy scales of the Hamiltonian for three different 
qubit systems, namely quantum dots, quantum shells, 
and Cooper-pair boxes, and we have calculated the mi- 
croscopic coupling constants J± and J z of the direct qubit 
interaction. 
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APPENDIX A: STRONGLY ASYMMETRIC 
QUBIT-FIELD COUPLING 

In order to clarify the physical processes that corre- 
spond to H ni i X in Eq. I|22|l . we consider the case of two 
qubits with identical level spacing, Q± — Q2 — ^> but 
with strongly asymmetric qubit-cavity mode coupling, 
gi = g 7^ 0, gi = 0. As will be discussed in more detail 
below, -f/mix then describes the competition between the 
direct coupling and the cavity-induced Stark and Lamb 
shift. For strongly asymmetric coupling, only qubit 1 



will experience a Stark and Lamb shift and hence will be 
shifted out of resonance with qubit 2. In particular, if 
the Stark shift of qubit 1 is large compared to J±, the 
qubits will become effectively decoupled. 

In order to quantify our argument, we evaluate the 
terms proportional to |V's,a)(V ; s,a| in Eqs. H2U|) and (|22|) in 
the limit of weak direct inter-qubit coupling, | Jj_\ <C \oj — 
(f2 ± J z /2)\. Then the denominators in the expressions 
for 7^ [Eq. I|17fl ] can be replaced by E S:(l — E- — > w — 
(il - J z /2) and by E+ - E s>a ->• u> - (Q + J z /2). Then, 



-hg* 



(qta + 1) ir+xr+i 

lo - (Q + J z /2) 



(Al) 



-hg' 



-Kg 



a^\T-)(T_\ 

LU-(fl- J z /2) 

j ataiitXitl 
u- (n + j z /2)' 



- fog ' 



(atq + i)| u)(tl 

w-(TL- J z /2) 



which represents the Stark and Lamb shift of qubit 1 . In 
particular, the last two terms in Eq. I|A1(I represent the 
relative energy shift of | ]{) and | {]) induced by the 
cavity. As discussed in Sec. IIIII this energy shift can be 
utilized to tune the qubits into and out of resonance by 
controlling the photon number in the cavity mode. 



APPENDIX B: SCHRIEFFER- WOLFF 
TRANSFORMATION FOR GENERAL 
HAMILTONIAN WITHOUT RWA 

For the coupling of the two-qubit system to the cavity 
via linearly polarized transitions, the full Hamiltonian 
reads in general 

flfiin = h J2 [ Qi§i >* + & ( a + at ) ( S t + 87)] ( B1 ) 

i=l,2 

+huja)a H (st§2 + H^t ) + hJ z S\ }Z S2, z - 

The energy non-conserving terms are usually neglected 
because of their smallness, leading to the Hamiltonian 
Eq. (Q. In this Appendix we analyze the additional 
contribution of the energy non-conserving terms to the 
coupled-qubit dynamics. 

As shown in Section lll CI the eigenstates of the coupled 
qubits are given by Eq. IjlUI) for vanishing qubit-cavity 
coupling. In this basis, 

i— ± i—s,a 

+/i(a + at)[ 7l |^)<T_|+7 2 |^ a )(T_| 

+ 73 |T + )<^| +74|T + )(^ a | + H.c] , (B2) 

with the coupling constants 7$ given by Eqs. I|13fl an d 
(|14fl . To first order in 7i , the Hamiltonian Hf u \i is diag- 
onalized by the Schrieffer- Wolff transformation 



lull 



iHtxe-^-ty+W-t+K (B3) 
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where B is defined by Eq. I|16|l and 

B+ = 7+|^)<r_|+7+|V>„)<T_| (B4) 
+7 3 + |T + )(^|+7 4 + |T 1 + )^a|, 

with the coefficients 



it = 




7i 




U) - 




E_ 


it = 




72 




LU - 


'rE a - 


E- 


it = 




73 




10 - 


IE+- 


■E. 


it = 




74 




to - 


h£+- 


E a 



(B5) 



Clearly, \^ff\ < for all i and u> > 0. Yet, 7+ ~ 7, for 
to < £7„ - B m , and |7+| ~ pftl for lu ^> E n - E m , where 
n m and n, m G {s, a, +, — }. 

Instead of Eq. (|TB|l , we obtain to second order in 7, the 
effective Hamiltonian, 



fjl tt 1 fVfull 1 rrfull , fVfull (TIA\ 

^full - H + S,L + + -^mix- (tfO) 



Here, Ho is given by Eq. Ijl9(l . The term for the Stark- 
and Lamb-shift is now calculated to be 



H s"l = J (« 2 + a t2 + 2a f a + 2) (B7) 

x {[7i(7i+7i + )-73(73 + 7 3 f )] lV>sX^I + [72(72 + 7^) - 74(74 +74~)] \^a)(ipa\ 
+ [73(73 + it) + 74(74 + 7 4 + )] IT+XT+I - [71(71 + 7i + ) + 72(72 + it)] |r_)<T_|} 

+H [73731^X^1 + 74 74 1 1pa)(lpa I + (7l7l +I2I2) \T-)(T-\] . 



In contrast to the result using the RWA, the additional 
terms containing 7^ modify the Stark shift, see Fig. [31 
and also the Lamb shift. In Fig. [3] we have assumed a 



photon Fock state in the cavity. Further, the two-photon 
transition term is given by 



#I- P h = \ [(737i " 7173 + 7472 - 7274) (a 2 |T+)(T_| + a^\T^)(T + \) (B8) 
+ (itli - 7i + 73 + ltl2 - itli) (a^ 2 \T + ){T_\ + a 2 \T_)(T + \) 
+ (7271 - 7374 + it 12 ~ it I3) (a t2 |^)(-0a| + a 2 \tp a )(tp s \) 
+ (11I2 - 7473 + 72~7i - 73" 74) (a 2 \tp s ){i> a \ + a) 2 \ip a ){ii) s \)] , 



which consists of i?2-ph and additional energy non- between \ip s ) and \ip a )- Similarly, 
conserving transitions between |T + ) and |T_) as well as 



^mt = \{ [7i(72 + it) + 72(7i + 7i + ) " 73(74 + it) ~ 74(73 + 7 3 + )] (« f « + + IVV.XV.I) (B9) 

+ (71 72" ~ 727l f + 7374 + 7473) IV'sXV'ol + (727l" - 7172 + 7374 + 7473) |^a)<^ S | 

+ [7i(73 + it) + 72(74 + 7 4 + ) - 73(7i + 7i + ) - 74(72 + it)] (« f « + 1)(|T+)<T_| + |T_)(T+|) 
+ (7371 + 7472 ~ 7i73 f ~ 7274 ) l T +X T -l + (717! + 737i + 727^ + HI2) \T-){T+\} 

I 



contains -ff m i x and new terms that couple |T+) and |T_) in addition to |-0 S ) and |^ Q ). In contrast to Eq. i|22[l . 
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the above term induces mixing of \T + ) and |T_) even for identical qubits. 
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